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1 Introduction 

As it was shown at first in pP, generalized 
cosmological Friedmann equations (GCFE) de- 
duced in the frame of gauge theories of grav- 
ity (GTG) |2] besides cosmological solutions 
lead to solutions for some hypothetical ob- 
jects — so-called superdense gravitating sys- 
tems (SGS). SGS have extremely high energy 
density and their dynamics is essentially none- 
insteinian. In considered approximation of ho- 
mogeneous isotropic space SGS oscillate be- 
tween minimum and maximum values of the 
scale factor. Later solutions for such systems 
in the case of various equations of state for 
gravitating matter were discussed in [3-5]. Be- 
low SGS are investigated for systems including 
scalar fields. 



2 Generalized cosmological 
Friedmann equations 

GCFE for systems including scalar field eft min- 
imally coupled with gravitation and radiation 



have the following form [7] 1 . 



R 2 



H 



Z 2 + \H 



hi 1 - 2/3(2F + j> 2 ) -3pV'4>} 



8tt 
3M 2 



1 - 2f3(2V + 



pUv 



+ H 

+ \2f3HcpV' 
-3(3 



Z 



8tt 
3M 2 



Pr + =<r + v 
z 

4(3(V - 40 2 ) 
1 - 2p(2V + ft 
-V' 2 )Z + 

V-tf-p r - 



(1) 

z, 



18/3 2 4 } 



(2) 



{v"4> 2 



2y,2 



^(3{W-tf) 2 



Z, 



where R is the scale factor, k = +1,0,-1 for 
closed, flat and open models respectively, H is 
the Hubble parameter, p r is radiation energy 
density, V(<j)) is a scalar field potential, V' = 

*v . " % Mp „ Planckain mass , z _ 



V" 



1 — /3(4V — 4> 2 ). (The system of units with h = 
c = 1 is used.) Eqs. (1)— (2) include indefinite 



1 One supposes the interaction between scalar field 
and radiation is negligibly small. 
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parameter (3 with inverse dimension of energy 
density, the value of \f3\~ 1 determines the scale 
of extremely high energy densities. 

Besides Eqs. (l)-(2) gravitational equa- 
tions of GTG lead to the following relation for 
torsion function S and nonmetricity function 

Q HE1 
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The conservation law in usual form follows 
from Eqs. (l)-(2), as result we obtain the 
equation for scalar field 



■ + 3H<j> = -V 



(4) 



and the integral for radiation p T R 4 = const. If 
the condition 



\f3(4V-<p 2 )\ «1 



(5) 



is valid, solutions of GCFE coincide practically 
with that of general relativity (GR) and the 
torsion and nonmetricity functions are negligi- 
bly small. 

As it was shown in [Jj, the GCFE allow to 
build regular inflationary cosmological models 
with dominating ultrarelativistic matter at a 
bounce, if the scale of extremely high energy 
densities is much less than the Planckian one 
\f3\~ 1 < 1M£ (P < 0). Numerical solutions ob- 
tained in |7J |S] for closed and flat models show 
that after inflation < lM p ) during cer- 
tain time intervals the dynamics of inflation- 
ary models has oscillating character, namely 
scalar field and the Hubble parameter oscillate 
near values zero. Such situation is typical for 
SGS. This means, if \cj)\ < lM p and the rela- 
tion (5) is not fulfilled, by choosing some initial 
conditions we can obtain numerical solutions of 
GCFE for SGS including scalar fields. We have 
to take into account that by given initial condi- 
tions for magnitudes (</>, (j),R, p r ) there are two 



different solutions corresponding to two values 
of the Hubble parameter, which according to 
(1) are 
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Solutions of GCFE corresponding to H + and 
H- were called in [Hj as i/ + -solutions and H-- 
solutions respectively. 

Initial conditions will be chosen at the 
moment t = corresponding to extremum of 
the scale factor R: R(0) = Ro and H + (0) = 
or -H-(O) = 0. Then according to (1) initial 
values of Rq, p r o, 4>q and 4>q satisfy the follow- 
ing relation 
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Eq. (6) determines on the plane P with the 
axes (0, 4>) curves (ifo-curves) , in points of 
which the Hubble parameter vanishes for H+- 
solutions or //--solutions. From (2) the time 
derivative of the Hubble parameter Hq at the 
moment t = is 
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where Zq = 1 — P(AVq — iy j 

In order to analyze solutions properties 
near the origin of coordinates on the plane P 
we will examine HQ-curves defined by (6) for 
some simple scalar field potentials. 
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3 Extremum i/g-curves for 
simplest scalar field poten- 
tials 

According to (1) in the case of flat and closed 
models we have Z > (/? < 0) and admissible 
values of scalar field eft and derivative 4> on the 
plane P are limited by two bounds L± defined 
by equation 

^ = ±{w + \p\- 1 y. (8) 

.f/o-curves B\ and B2 for flat models (k = 0) 
are situated on the plane P between bounds 
L + and L_ having with them common points 
(points K\ and K2 in Fig. 1) on the axis 
(j) (if V'(0) = 0). Each of HQ-curves con- 
tains two parts — extremum curve for H + - 
solutions (Ho + -curves) and extremum curve 
for ff_-solutions (i?o-- curves ) |Ej- The family 
of HQ-curves for closed models is situated on 
the plane P between curves B\ and i?2- The 
family of -ffo _curves f° r open models is situated 
on the plane P between bounds L± and corre- 
sponding f/o- curves for flat models. 2 

For given scalar field potential V{(f>) the 
behaviour of iio- curves depends on parame- 
ters of potential V, on parameter (3 and on 
the model type (k = 0, ±1). Below Hq are 
analyzed for scalar field potential V = ^m 2 ^ 2 
(rn = 10 _6 Mp) in the case of p r = 0. For 
flat models -£fo- curves are given in Fig. 1 (/? = 
— 1.6 • 10 12 ). The derivative Hq is positive in 
points of -ffo-curves if the value of \/3\ < \f3±\ ~ 
1.8 • 10 12 . This means that .ffo- curves i n con- 
sidered case are bounce curves discussed in jSj. 
By increasing of the value of \(3\ the points 
K\ and K2 approach to the origin of coordi- 
nates and the derivative Hq becomes negative 
in the neighbourhood of points K\ and K2 (see 
Fig. 2). In points of -£/o- curves with Hq < the 

2 In the case of open models k = — 1 Eq. (6) allows 
solutions if Z < 0. These solutions are not discussed in 
present paper. 
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Figure 1: i/o-curves of flat models at (3 = —1.6- 
10 12 
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Figure 2: i^o-curves of flat models at = —1.0- 
10 14 
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Figure 3: iio-curves of closed models at f3 
-1.0 • 10 14 



evolution of system corresponds to transition 
from expansion to compression. The existence 
of regions on iio-curves with Hq < is nec- 
essary condition for appearance of oscillating 
solutions for SGS. The behaviour of iio-curves 
for open models is similar to that for flat mod- 
els. The curves B\ and B2 are limiting in- 
curves for open models under Rq — > 00. By de- 
creasing of the value of Rq corresponding in- 
curves of open models approach to bounds L± . 
All iio-curves of open models contain points 
K\ and K2, where H + = Like to flat 

models oscillating SGS-solutions of open type 
can appear by sufficiently large values of \(3\. 

In the case of closed models by certain 
value of (3 there are three kinds of incurves 
(Fig. 3). For very large values of Ro > Rqi 
iio-curves are closed curves with the center in 
origin of coordinates (curve (a) in Fig. 3) , they 
correspond to iio-curves of GR. The value of 
i?oi depends on \f3\, for example i?oi ~ 7.14 • 
1.0 • 10 14 M- 4 . In the case 



of i?o < R01 we have i^o-curves (6) and (c) 
presented in Fig. 3. The Hubble parameter in 
points of curves (b) is positive (bounce-curves), 
and it is negative in points of curves (c) . So, in 
points of iio-curves (</>o 0, 4>o = 0) the values 
of i?o and Vq are connected by usual relation of 

Vq, and the derivative iio = 
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The presence on iio-curves of regions with dif- 
ferent signs of H + (H-) leads to appearance 
of oscillating solutions. Note, that in points of 
iio-curves lying on coordinates axes we have 
H+ = 

The analysis shows that the behaviour of 
iio in the case of potential V2 = \X4> 4 (A = 
const) qualitatively is like to that of considered 
above potential V\ = \m 2 4> 2 . 

4 Numerical analysis of solu- 
tions for scalar fields SGS 

By using various scalar field potentials apply- 
ing in inflationary cosmology (in particular, 
V\ = ^m 2 (j) 2 , V2 = jA0 4 ), numerical analysis 
of SGS solutions of Eqs. (1), (2), (4) for flat, 
open and closed type was carried out, solutions 
properties in dependence on value of parameter 
P and initial conditions were investigated. As 
result the following conclusions were obtained: 

1. SGS-solutions of flat and open type have 
non-stationary character and are unsta- 
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ble. Particular SGS-solution of flat type 
is given in Fig. 4. 3 

SGS-solutions of closed type can have 
quasi-stationary as well as non-stationary 
regime in dependence on the value of 
(3 {w - $j . The value of 4> 2 has to 
be less than 4Vo at several orders to con- 
serve the stable character of discussed so- 
lutions. In connection with this the anal- 
ysis of SGS-solutions will be given below 
in dependence on the value of x = 4|/3Vq|. 



3. Quasi-stationary regime for SGS-solutions 
of closed type takes place, if the value 
x ~ 1 and changes in certain interval 
a^min < x < x max depending on the value 
of (3 and the potential form. So, in the 
case of potential V\ we have 1.5 < x < 4.5 
at |/3| = 1.0- 10 14 . By increasing of \fi\ the 
values of x m i n and x max decrease. For ex- 
ample, in the case of V\ and \(3\ = 1.0- 10 20 
we have 0.12 < x < 2.0. The frequency of 
(^-oscillations and also frequency and am- 
plitude of .ff-oscillations do not depend 
practically on the value of (3 for given 
form of scalar field potential in the case of 
quasi-stationary regime. The realization 
of quasi-stationary regime means, that the 
sign of term 3H(j) in Eq. (4) changes dur- 
ing (^-oscillations period, during one half 
period this term plays the role of damp- 
ing force and during the second half pe- 
riod it plays the role of accelerating force. 
As result, the frequency of if-oscillations 
is two times greater than the frequency 
of ^-oscillations. In accordance with (3) 
torsion (nonmetricity) function oscillates 
with the same frequency as the Hubble pa- 
rameter 4 . The evolution of magnitudes 4>, 



H, R and the energy density of scalar field 
Pcj> = ^ 2 + V(4>) is presented in Fig. 5. 
Note, that scalar field energy of closed 
SGS changes with the time also. 

4. By decreasing of the value of x ra \ a (x < 
%mm) oscillating SGS-solutions of closed 
type become unstable. 

5. If x S> 1 oscillating SGS-solutions of 
closed type are essentially non-stationary 
for limited time intervals (Fig. 6). In the 
case of large time scale oscillations have 
quasi-modulated character. 

6. In accordance with the structure of GCFE 
the presence of radiation has essential in- 
fluence on dynamics of SGS in the case 
of quasi-stationary regime x ~ 1 destroy- 
ing this regime, and such influence is es- 
sentially less in the case of non-stationary 
regime x 1. 

In conclusion note, that oscillating classical so- 
lutions for scalar fields are unstable with re- 
spect to quantum processes of transformation 
of scalar fields into elementary particles and 
hence such solutions can have physical sense 
only for limited time intervals. 
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4 Quasi-stationary oscillating solutions for closed 
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singular cosmological model. From our point of view 
satisfactory cosmological model can not be built by us- 
ing solutions for scalar fields SGS. 
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Figure 4: SGS-solution of flat type (/? = -1.0 • 
10 20 , <f) = 5.25 • 1(T 5 and O = 7.09 • 10~ 12 ). 
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Figure 5: SGS-solution of closed type in quasi- 
stationary regime (/? = — 1.0 • 10 14 , 0o = 0.1 
and 4>q = 0). 
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Figure 6: SGS-solution of closed type in non- 
stationary regime (/? = — 1.0 • 10 16 , <p$ = 0.1 • 
10~ 5 and 0o = 0.) 



7 



